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We present an analytical approach to the theory of nonlinear propagation of femtosecond opti¬ 
cal pulses with broad-band spectrum in gases. The vector character of the nonlinear third-order 
polarization is investigated in details, taking into account the carrier to envelope phase. The corre¬ 
sponding system of vector amplitude equations is written by using left-hand and right-hand circular 
components of the electrical field. We found that this system nonlinear equations admits 3 D + 1 
vector soliton solution with Lorentz shape. The solution presents relatively stable propagation and 
rotation with GHz frequency of the vector of the electrical field in plane, orthogonal to the direction 
of propagation. The evolution of the intensity profile demonstrate weak self-compression and week 
spherical wave in the first milliseconds of propagation. 

PACS numbers: 42.65.Sf, 52.38.Hb 


INTRODUCTION 


When a femtosecond laser pulse with power above 
the critical for self-focusing propagates in air, a num¬ 
ber of new physical effects have been observed, such as 
long-range self-channeling |l, 4], coherent and incoherent 
GHz and THz emission [3, 4], asymmetric pulse shap¬ 
ing, super-broad spectra [S, (6], sell-compression [§j], po¬ 
larization rotation [j] and others. A remarkable effect 


is also that some of the light pulses propagate over dis¬ 
tances of several kilometers in vertical traces, preserving 
their spectrum and shapes 0. In one typical experi¬ 
ment in the near zone up to 1 — 3 m from the source, 
when the pulse’s intensity exceeds I > 10 12 Wycm 2 , ini¬ 
tial self-focusing and self-compressing starts, which leads 
to enlarging the k z spectrum to broad band asymmetric 
one A k z ~ k$. The process increases the core inten¬ 
sity up to 10 14 IU/cm 2 , where a short non-homogenous 
plasma column in the nonlinear focus is observed. Usu¬ 
ally the standard model describing the propagation in the 
near zone is a scalar spatio-temporal nonlinear paraxial 
equation including in addition terms with plasma ion¬ 
ization, higher order Kerr terms, multiphoton ionization 
and others 1,0- The basic model works properly in 
the near zone because of the fact, that paraxial approx¬ 
imation is valid for pulses with narrow-band spectrum 
A k z « k 0 . In the far-away zone plasma generation and 
higher-order Kerr terms are also included as necessary for 
the balance between the self-focussing and plasma defo¬ 
cussing and for obtaining long range self-channeling in 
gases. However, the above explanation of filamentation 
is difficult to be applied in the far-away zone. As reviewed 
Mi plasma density at long distances from the 


source is to small to prevent self-focusing. There are basi¬ 
cally three main characteristics which remain the same at 
these distances - the broad-band spectrum, coherent GHz 
generation and the width of the core, while the inten¬ 
sity in a stable filament drops to a value of 10 12 IU/cm 2 . 


The higher-order Kerr terms for pulses with intensities 
of order of / ~ 10 12 IU/cm 2 are also too small to pre¬ 
vent self-focussing. The experiments, where observation 
of long-range self-channeling without ionization was re¬ 
alized II413I, show the needing to change the role of 
the plasma in the laser filamentation. In addition, there 
are difficulties with the physical interpretation of the co¬ 
herent GHz radiation as a result of plasma generation. 
The light pulse near the nonlinear focus emits incoherent 
and non-homogenous plasma, while the coherent GHz ra¬ 
diation requires homogenous plasma with fixed electron 
density of the order of 10 15 cm -3 . Only homogenous 
plasma can generate coherent GHz emission, but such 
kind of plasma absent in the process of filamentation. In 
the experiments the following basic characteristic of the 
long-range filaments in air are detected: 

1. Broad-band spectrum (A k z ~ /co). 

2. Intensity a little above of the value for self-focusing 
{I — heif-foc — 10 12 W/cm 2 ). 

3. Absence of plasma at long distances. 

4. Asymmetric relatively stable (Lorentz) spectral and 
longitudinal forms. 

5. Coherent GHz generation. 

Recently in was developed scalar ionization-free 
non-paraxial nonlinear model which gives the above char¬ 
acteristics of the stable filament. The analytical and the 
numerical results describes correctly the linear and non¬ 
linear evolution of narrow-band and broad-band laser 
pulses and in addition it was found, that the equation 
admits exact Lorentz-type soliton solutions in approxi¬ 
mation of neglecting the GHz oscillation. Nevertheless 
this theory encountered significant difficulties. The main 
problems are: 

1. Peak instability of the soliton solution under small 
initial perturbations. 

2. The soliton solution is obtained after neglecting the 
GHz oscillation. 

3. The soliton solution admits one free parameter. 
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4. There are problems with the conservation law of the 
nonlinear operator, when we use the GHz oscillation. 

To solve the above problems we propose in this paper 
a nonlinear vector model. 


components are orthogonal and complex conjugated. In 
the case of Maker and Terhune polarization © we obtain 
for A x and A y components the following expressions: 


NONLINEAR POLARIZATION 

The self-action process broadens the pulse spectrum - 
starting with narrow-band pulse, the stable filament be¬ 
comes broad-band far way from the source. From other 
hand, the standard filamentation model based on plasma 
generation and multi-photon processes includes also non¬ 
linear polarization of the kind 


P nl = n 2 


(e-e*^e + ^(e- Ej E* 


(1) 


where n 2 is the nonlinear refractive index of the isotropic 
media. The polarization © was proposed by Maker and 
Terhune in 1965 14] . If the electrical field contains one 
linear or circular component, the polarization (|TJ) de¬ 
scribes only the self-action effect, while in the case of 
two-component electrical filed E = (E X) E y , 0) additional 
terms appear, presenting cross-modulation and degener¬ 
ate four-photon parametric mixing. The self-action pro¬ 
cess broadens the pulse spectrum - starting from narrow- 
band pulse, the stable filament becomes broad-band far 
way from the source. Later in 15-17] it is shown, that the 
evolution of broad-band pulses like filaments can not be 
described correctly by nonlinear polarization of the kind 
([T]). It is more correct to use the generalized nonlinear 
operator 


P nl =n 2 (E- E) E, (2) 

which includes additional processes associated with third 
harmonic generation (THG). The more precise analysis, 
presented in the paper, demonstrates that the polariza¬ 
tion of kind © is not applicable to a scalar model, be¬ 
cause the corresponding Manley-Rowe (MR) conserva¬ 
tion laws are not satisfied. That is why we substitute 
into the nonlinear operators © and © two-component 
electrical vector field at one carrying frequency in the 
form: 


p? = >2 


= S" 2 


I-4.I 2 + jl4l 2 ) -4. - l^Al 


x exp [iko(z — v p ht )] + c.c. 


(4) 


\A y \ 2 + gl^l 2 ') Ay ~ 3 A*Al 


x exp [iko(z — v p ht )] + c.c. 


The nonrestricted nonlinear polarization © generates 
the following components 


P.f = n 2 


(A 2 - A 2 y ) A x exp [2ik 0 (z - v ph t)} 


+ l\A x \ 2 + -\A y \ 2 \A x --A* x A 


exp [■ ik 0 {z - v ph t)] + c.c. 

(5) 


pnl _ ^2 

y ~ i 


( A 2 - A 2 ) Ay exp [2ik 0 (z - v p ht)\ 


1 


+ l\A y \ 2 + =-\A x \ 2 )A y -^A*A 2 


exp [ ik 0 (z - v ph t)] + c.c., 


where h 2 = | n 2 . Comparing © and ©, it is clearly 
seen that the operator © n 2 (^E ■ E'j E generalizes the 
case of Marker and Terhune’s operator ©, but includes 
also additional terms associated with THG. 

The generalized nonlinear polarization © is quite sim¬ 
ple in the terms of left-hand and right-hand circularly 
components fl9l. [20| 


E + (x,y, z,t) = A + (x,y 1 z,t)ex p [ik 0 (z - v ph t)] 


( 6 ) 


E- (x, y,z,t) = A_ (x, y,z,t) exp [~ik 0 (z - v ph t)\ 


where 


- _ {A x exp [ik 0 (z - v ph t)] + c.c.) _ 

lA — -—- X - 


(Ay exp [ik 0 (z - v p ht)] - c.c.) 


2 i 


(3) 


A+ — (A x + iA y ) / \[2 
A- = (A x — iAy)/y/ 2, 


(7) 


where A x = A x (x, y , z, t),A y = A y (x, y, z, t ) are the am¬ 
plitude functions, ko is carrying wave number, uj o is the 
carrying frequency of the laser source and v p h is the phase 
velocity. Here we take in mind the fact that the x and y 


are the left-hand A + and right-hand A_ circular com¬ 
ponents of the amplitude fields. Using the technique of 
calculations, presented in fl9| . we obtain that the nonlin¬ 
ear polarization © generates the following components 
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P+ = n 2 ( A 2 + A _) exp [ik 0 (z - v ph t )] 

P- = n 2 (A 2 _A + ) exp [-ik 0 (z - , 

where n 2 is the nonlinear refractive index. 


( 8 ) 


BASIC SYSTEM OF EQUATIONS 

The dynamics of narrow-band laser pulses can be ac¬ 
curately described in the frame of paraxial optics. The 
filamentation experiments demonstrate a typical pulse 
evolution: the initial laser pulse (to > 50/s) possesses a 
relatively narrow-band spectrum (A k z -C fco) and during 
the process the initial self-focusing and self-compression 
broadens the spectrum significantly. The broad-band 
spectrum (A k z ~ fco) is one of the basic characteristics of 
the stable filament. The evolution of the filament can not 
be further described in the frame of the nonlinear parax¬ 
ial optics because the paraxial optics works correctly for 
narrow-band laser pulses only. 

Recently, the evolution of broad-band laser pulses was 
described correctly within different non-paraxial models 
such as UPPE 15] or non-paraxial envelope equations 
[lij ]. Another standard restriction in the filamentation 
theory is the use of one-component scalar approximation 
of the electrical field E. This approximation, though, is 
in contradiction with recent experimental results, where 
rotation of the polarization vector is observed [7|. For 
this reason in the present paper we use the non-paraxial 
vector model up to second order of dispersion, in which 
the nonlinear effects are described by the nonlinear polar¬ 
ization components ([ 8 ]). The system non-paraxial equa¬ 
tions of the amplitude functions A + and A_ has the form 


— 2* fco 



1 dA + \ 

Vgr dt J 


2iko 


dA_ 1 dA_ X 

dz Vgr dt J 


= AA 


1 + /3 d 2 A + 

' W 

+fco?7, 2 A]j_A_ 


(9) 


= AA_ - 


1 + /3 d 2 A_ 

~^f r W 


A fc o n 2 A 2 _ A . 


where A is 3 D - (x, y , z) Laplace operator, v gr is the 
group velocity, /3 = kov 2 r k" and fc" is the group velocity 
dispersion. In the process of deriving of the system of 
Eqs. ©I is used in mind, that the left-hand and right- 
hand components are complex conjugated fields. 

This model describes the ionization-free filamentation 
regime and we will investigate the case, when the pulse 
intensities are slightly above to the critical one for self- 
focusing. The dispersion number in air is too small 


/3 = kov 2 r k" ~ 2.1 x 10 5 and we will use Eqs. © 
in approximation of the first order of dispersion. 


— 2ifco 


rM + | 1 

dz Vgr dt J 


= AA+ 


2iko 



1 dA. \ 

- ^r) =AA 

Vgr dt J 


1 d 2 A+ 
v 2 gr dt 2 
+7 A 2 + A_ 

( 10 ) 

1 d 2 A_ 

Vgr dt 2 

+7 A 2 _A+, 


where 7 = n 2 fcQAQ is nonlinear coefficient. 


VECTOR SOLITON SOLUTION AND VECTOR 
ROTATION. THE FILAMENT AS ROGUE WAVE 

The nonlinear system of equation m admits exact 
solitary vector solution, when 7 = 2 and the spectral 
width of the pulses reach the value Afc z a: fco 

2 

A+(x, y, z, t ) = -—37 exp [-*Afc z (z - v gr t)] 

1 + r- 

( 11 ) 

2 

A- (x , y , z, t) = 1 + _ 2 exp [iAk z (z - v gr t )], 


where f = y x 2 + y 2 + (z - ik cep ) 2 - v 2 r (t - ik ce p/vgr) 2 
and k cep is determinate below. The solution of the corre¬ 
sponding vector electrical field can be written after multi¬ 
plying the amplitude functions m with the main phase 

m 


E + (x , y, z, t ) = exp [-*Afc 2 ( v ph - v gr ) t] 

( 12 ) 

2 

£_ (x, y,z,t) = ^ _ 2 exp [iAk z (v ph - v gr ) t\. 

Lets turn from the left-hand and right-hand circular com¬ 
ponents 0 to the standard Cartesian coordinates 

E x = (. E+ + E-)/\/2, E y = (. E+ - E_)/(iV 2). (13) 

The solution (fT21) written in Cartesian coordinates has 
the form 


y,Z,t) — COS [Afc z ( Vph Vgr^j t\ 

1 + r z 

(14) 

2 

E v (x, y, z, t) = 1+ ~ 2 sin [Afc z (v ph - v gr ) t]. 
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a) 


b) 


c) 




z=0 


z=20 cm 


z=10m 


FIG. 1. Evolution of the intensity profile I = |Ae| 2 + |^4^| 2 of the solution (11411 . Side (x, z) projection is plotted. On Fig.la) 
is the initial profile, while on Fig.lb) the evolution of the intensity profile on distance z=20 cm is presented. It is clearly seen 
the generation of weak spherical wave. At large distances the solution is relatively stable, as it can be seen from the Fig.lc). 
One more detailed picture of the propagation in the far-away zone demonstrate a weak self-compression of the pulse. 


The 3 D + 1 Lorentz type solution (ED. presented in 
Cartesian coordinates, gives oscillation of the electrical 
vector E = (E x ,E y , 0) in [x,y) plane. It can be seen 
directly that the frequency of oscillation is equal to the 
carrier to envelope ones uj cep = ko ( v p h — v gr ). The cor¬ 
responding longitudinal spatial carrier to envelope wave 
number is k cep = u> cep /v gr . In standard atmosphere this 
complies to a rotation of the electrical field vector with 
period in z direction X cep = 10 — 20 cm. Other impor¬ 
tant characteristic of the vector solution ED is that it 
is soliton in 3 + 1 dimensional complex Minkowski space 
and his real dynamics in 3D space is more close to weak 
Rogue wave with small self-compressing in the direction 
of propagation. The evolution of the intensity profile is 
plotted on Fig.l and demonstrate a generation of week 
spherical wave (Fig 1. b) in the first milliseconds of the 
propagation distance. At large distances, as it can be 
seen from the Fig. lc), the propagation is relatively sta¬ 
ble with preserving the amplitude maximum. One more 
detailed picture of the propagation in the far-away zone 
demonstrate a weak self-compression of the pulse. 


CONCLUSIONS 

In this investigation we point that in the femtosecond 
region it is impossible to reduce the nonlinearity of third 
order to the scalar Kerr type only (proportional to the 
intensity). Therefore, we study more precisely the vec¬ 
tor character of the nonlinear third order polarization, 
taking into account the carrier to envelope phase. The 
corresponding amplitude equations are written in terms 
of left-hand and right-hand circularly polarizing compo¬ 


nents. The vector system ED admits exact (3+l)D soli¬ 
ton solutions El of Lorentz type. Our soliton solution 
is obtained for broad-band pulses A k z ~ kg and leads 
to the conclusion that the soliton appears as a balance 
between parabolic divergent type wave diffraction and 
parabolic convergent nonlinearity. The solution gives also 
a rotation of the vector of the electrical field with pe¬ 
riod equal to the period of carrier to envelope wavelength 
A cep = 10 — 20 cm, in depends from the air pressure and 
temperature. Other important result is the generation of 
weak spherical wave in the first milliseconds of propaga¬ 
tion and weak self compression in the far-away zone. 


* lubomirkovach@yahoo.com 
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